Image  Cover  Sheet 


CLASSIFICATION 


UNCLASSIFIED 


SYSTEM  NUMBER 


503453 


TITLE 


HYDRODYNAMIC  INTERACTION  FORCES  BETWEEN  TWO  SHIPS  IN  WAVES 


System  Number 
Patron  Number 
Requester : 


Notes : 


DSIS  Use  only: 

Deliver  to; 


HYDRODYNAMIC  INTERACTION  FORCES  BETWEEN  TWO  SHIPS  IN 
WAVES 


Y.  J.  He,  Centre  for  Marine  Vessel  Design  and  Research,  Faculty  of  Engineering,  Dalhousie  University, 
(formerly  Technical  University  of  Nova  Scotia),  Halifax,  Nova  Scotia,  B3J  2X4 

K.  McTaggart,  Defence  Research  Establishment  Atlantic,  RO.  Box  1012,  Dartmouth,  Nova  Scotia,  B2Y 
3Z7 

C*  C*  Hsiung,  Centre  for  Marine  Vessel  Design  and  Research,  Faculty  of  Engineering,  Dalhousie 
University,  (formerly  Technical  University  of  Nova  Scotia),  Halifax,  Nova  Scotia,  B3J  2X4 


Abstract 

The  hydrodynamic  interaction  forces  in  waves  are  computed  for  two  ships  advancing  in  par¬ 
allel  at  the  same  speed .  The  forces  considered  are  due  to  added  mass,  wave  radiation  damp¬ 
ing,  and  wave  excitation  including  diffraction  effects .  The  coupled  hydrodynamic  forces  are 
computed  by  the  three-dimensional  panel  method  using  the  zero  speed  Green  function,  with 
extensions  to  include  forward  speed  effects.  The  theoretical  formulation  correctly  satisfies 
boundary  conditions  on  both  ship  hulls  simultaneously .  The  wave  diffraction  problem  is 
solved  for  the  presence  of  both  ships  in  a  stationary  state.  The  wave  radiation  problem  is 
solved  by  assuming  one  ship  is  in  motion  while  the  other  is  at  rest.  Examples  are  given  for 
two  stationary  vertical  circular  cylinders  and  for  a  small  vessel  advancing  in  close  proximity 
to  a  large  supply  vessel. 


NOMENCLATURE 

a  wave  amplitude 

Bi  beam  of  ship-b 

Ff  force  component  j  on  ship-a 

Fj*a  radiation  force  component  j  on  ship-a 

Fffa  wave  excitation  force  component  j  on  ship-a 

fja  complex  incident  force  j  on  ship-a 

fDa  complex  diffracted  force  j  on  ship-a 

k  wavenumber 

Lb  length  of  ship-b 

rn%  generalized  normal  at  speed  for  mode  k 
of  ship-a 

Uj  normal  on  ship-a  for  mode  j 
p  hydrodynamic  pressure  or  field  point 

q  source  point 

R  horizontal  distance  from  field  point  to  source 

r  cylinder  radius 

ri  distance  from  field  point  to  source 

7*2  distance  from  field  point  to  image  source 

fa  moment  arm  vector 

Sa  surface  of  ship-a 

T  draft 

U  ship  speed 

x,y,z  coordinate  system 

Xbi  yb  position  of  ship-b 

W  steady  flow  velocity  vector 

/?  wave  direction 


Q  complex  motion  amplitude  mode  k  of  ship-a 
A  wavelength 

A||  damping  component  j  of  ship-a  for 
motion  of  ship-b  in  mode  k 
Pjk  added  mass  component  j  of  ship-a  for 
motion  of  ship-b  in  mode  k 
v  radiation  wavenumber 

p  water  density 

dp  diffracted  source  strength  on  ship-a 

a%b  source  strength  on  ship-b  due  to 

mode  k  of  ship-a 

$  unsteady  wave  velocity  potential 

<I>r  total  velocity  potential 

<j)n  complex  diffracted  wave  potential 

0/  complex  incident  wave  potential 

<f>%  complex  potential  for  mode  k  of  ship-a 

<()r  complex  radiated  wave  potential 

<j>8  disturbed  steady  potential 

oj  wave  frequency 

u)e  wave  encounter  frequency 

A  a  displacement  of  cylinder-a 


INTRODUCTION 

When  two  ships  are  advancing  together  at 
speed  in  close  proximity  in  waves,  the  hydrody¬ 
namic  forces  on  each  ship  are  affected  by  shelter- 
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ing  or  wave-reflection  from  the  neighbouring  ship. 
This  paper  describes  research  that  was  initiated 
by  the  design  process  for  a  naval  supply  vessel, 
for  which  hydrodynamic  interaction  effects  with 
receiving  ships  can  be  significant.  The  resulting 
computer  program  can  predict  added  mass,  wave 
radiation  damping,  and  wave  excitation  forces  act¬ 
ing  on  each  vessel  in  the  presence  of  the  other. 

The  present  code  is  an  extension  of  a  panel 
code  by  Huang  and  Hsiung  [1]  that  was  used  by 
Howard,  Calisal  and  Barnard  [2]  for  computing  hy¬ 
drodynamic  forces  and  resulting  motions  of  a  single 
ship  in  waves.  The  theory  for  incorporating  inter¬ 
action  effects  between  two  bodies  has  been  devel¬ 
oped  by  Van  Oortmerssen  [3],  and  Kokkinowrachos 
[4].  While  codes  for  the  single  ship  case  utilize 
lateral  symmetry,  solution  of  the  non-symmetrical 
two  ship  case  using  the  current  code  requires  ap¬ 
proximately  four  times  as  many  panels  as  for  a 
single  ship.  Consequently,  the  computational  re¬ 
quirements  for  memory  and  CPU  time  are  of  the 
order  of  sixteen  times  greater  than  for  an  individ¬ 
ual  ship  case. 

As  a  preliminary  validation,  computed  forces 
are  compared  with  published  numerical  and  exper¬ 
imental  results  by  Kokkinowrachos  [4]  for  two  ver¬ 
tical  cylinders  floating  in  waves  with  zero  forward 
speed.  To  illustrate  the  influence  of  interaction  ef¬ 
fects  on  ship  operations,  predicted  hydrodynamic 
forces  are  presented  for  a  small  vessel  abeam  of  a 
large  supply  vessel. 


FORMULATION 

Figure  1  shows  a  plan  view  for  solution  of  hy¬ 
drodynamic  interactions  between  two  ships.  The 
coordinate  system  has  its  origin  in  the  calm  water 
surface  plane  at  midships  of  ship-a.  The  coordi¬ 
nates  Xb,yb  represent  the  location  of  midships  of 
ship-b.  The  incident  wave  direction  /?  follows  the 
convention  of  180  degrees  for  head  seas  and  90  de¬ 
grees  for  waves  from  starboard.  The  £  axis  has  a 
positive-up  direction. 


Velocity  Potentials 


Figure  1:  Plan  View  of  Ship  Interaction 

be  expressed  as: 

$r(x,y,z,t)  = 

~Ux  +  <j)a(x,y,z)  +  $(x,y,z,t)  (1) 

where  U  is  steady  forward  speed  of  both  ships,  <f)8 
is  the  disturbed  potential  due  to  the  presence  of 
the  ship  in  the  steady  incident  flow,  and  <3>  is  the 
total  oscillatory  potential  due  to  waves  and  ship 
motions.  The  term  —  Ux  represents  the  incident 
potential  due  to  ship  forward  speed.  The  unsteady 
velocity  potential  <I>  can  be  written  as: 

$( x,y,z,t )  =  Re{[<j)I{x,y,z)  +  (pD(x,y,z) 

+<pR{x,y,z)]e~i“’t]  (2) 

where  Re  denotes  the  real  part  of  the  term  in 
brackets,  and  <j)i,  4>d  and  <f)R  are  the  complex  val¬ 
ues  of  the  incident,  diffracted,  and  radiated  poten¬ 
tials,  and  coe  is  the  wave  encounter  frequency. 


Hydrodynamic  Forces 

The  hydrodynamic  forces  acting  on  ship-a  and 
ship-b  can  be  expressed  as: 


F?  = 

j  prijdS 

U  =  1, 2, 

— >6) 

on  Sa  (3) 

II 

f  prijdS 

U  =  1)2,  • 

on  Sb  (4) 

Sb 


where  Sa  and  Sb  are  the  wetted  hull  surfaces  of 
ships  a  and  b,  and  p  is  hydrodynamic  pressure. 
The  generalized  unit  normal  of  each  ship  can  be 
expressed  as  shown  here  for  ship-a. 


For  oscillatory  forces  acting  on  a  body  in 
waves,  viscous  effects  are  small  and  it  can  be  as¬ 
sumed  that  the  fluid  is  inviscid  and  incompressible 
and  that  the  flow  is  irrotational.  The  resultant  ve¬ 
locity  potential  4>y  in  the  vicinity  of  the  ships  can 


f  n j  for  j  =  1, 2, 3 

\  ra  x  for  j  —  4, 5, 6 


(5) 
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where  ra  is  the  moment  arm  vector  from  the  origin 
to  the  point  on  the  hull.  For  j  —  1,  2, 3,  n't  are  the 


components  of  the  normal  vector  pointing  into  the 
hull  of  ship-a. 

The  hydrodynamic  pressure  p  can  be  com¬ 
puted  from: 

P  =  -p(^  +  ^-V$)  (6) 

where  p  is  water  density.  The  steady  flow  velocity 
vector  W  is  given  by 

W  =  V(-*7x  +  0S)  (7) 

Substituting  Equations  (6)  and  (7)  into  Equa¬ 
tions  (3)  and  (4),  the  hydrodynamic  force  on  the 
two  ships  can  be  expressed  as: 


n 


-M 

Jsa  L 


d$_ 

dt 


+  V{-Ux  +  <j>s)  ■  V$ 


n]dS 


Fi  = 


~  f  P  Sr  +  V(-J7*  +  4>.)  •  V$  n]dS 
Jsh  lut  J 


(8) 

(9) 


By  substituting  the  various  components  of  Equa¬ 
tion  (2)  into  the  above  two  equations,  one  can  com¬ 
pute  the  incident  wave  force,  diffracted  wave  force, 
or  radiated  wave  forces  acting  on  each  ship. 


RADIATION  PROBLEM 


for  unit  amplitude  motions  of  ship-a  while  ship-b 
is  at  rest,  and  ££  and  4>bk  are  corresponding  values 
for  ship-b. 

The  radiation  velocity  potentials  <pk  and  4>\ 
must  satisfy  several  conditions.  The  solution  of  <f>% 
to  satisfy  these  conditions  will  be  presented  here, 
with  a  similar  solution  existing  for  (j)bk. 

For  potential  flow,  the  radiation  potential  <j)k 
must  satisfy  Laplace’s  equation  throughout  the 
fluid: 


V2<Pt  =  0  (11) 

The  following  boundary  condition  applies  on  the 
mean  free  surface: 

(9iL + u2ih + 2iuJeUik  ~  ^ 

=  0  on  z  ~  0  (12) 

The  body  normal  boundary  conditions  on  the  mov¬ 
ing  ship-a,  the  stationary  ship-b,  and  ocean  bot¬ 
tom  are: 


m 

dn 

=  K 

mk  a 

—  — ~  On  Sa 

l(Je 

(13) 

dn 

=  0 

on  Sb 

(14) 

m 

dz 

=  0 

on  z  —  —co 

(15) 

The  solution  of  hydrodynamic  forces  due  to 
ship  motions  is  called  the  radiation  problem.  For 
the  present  case  of  two  ships  in  close  proximity,  the 
radiation  solution  must  include  interaction  effects. 

Radiation  Potential 

For  two  ships  freely  floating  in  calm  water, 
the  radiation  wave  velocity  potentials  can  be  de¬ 
termined  by  satisfying  the  body  surface  condition 
of  two  cases:  1)  ship-a  is  in  motion  while  ship-b  is 
at  rest;  and  2)  ship-b  is  in  motion  while  ship-a  is 
at  rest. 


where  mk  is  the  generalized  body  normal  vector 
including  the  effect  of  ship  forward  speed,  as  de¬ 
scribed  in  Reference  1.  The  velocity  potential  <pk 
must  also  satisfy  the  condition  that  waves  radiate 
outward  from  the  body. 

The  solution  of  the  velocity  potential  <f)k  is 
facilitated  by  a  Green  function  that  satisfies  all 
boundary  conditions  except  the  body  boundary 
conditions  of  Equations  13  and  14.  These  bound¬ 
ary  conditions  can  be  satisfied  by  placing  sources  of 
strength  uka  on  ship-a  and  strength  akb  on  ship-b. 
The  radiation  potentials  due  to  motions  of  ship-a 
while  ship-b  is  stationary  can  then  be  expressed  as: 


The  complex  radiation  potential  from  ships  a 
and  b,  each  having  six  degrees  of  freedom,  can  be 
expressed  as: 

6 

4>r(*,  y , *, *)  =  £  [CM  +  CM]  (io) 

k- 1 

where  ££  is  the  complex  motion  amplitude  of  ship- 
a  for  mode  k ,  <f>k  is  the  complex  velocity  potential 


«(P)  =  /  <a(q)G(p;q)dS(q) 

+  [  ^b(q)G{p;q)dS(q)  (16) 
Jsh 

where  p  is  the  location  of  a  point  in  the  fluid 
domain,  G(p;q )  is  the  Green  function,  and  q  is 
the  location  of  a  point  on  the  hull  surface.  The 
Green  function  satisfying  the  continuity  equation 
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and  most  of  the  boundary  conditions  is: 


11  f°° 

G(p]  q)  = - 1 - b  2is  /  - JQ{uR)dk 

n  r2  Jo  k  —  v 

+i‘l-Kvev^*+z^  Jo{yR)  (17) 

where  xp,yp,zp  are  the  coordinates  of  point  p, 
xq<>  Vqi  zq  axe  the  coordinates  of  point  q ,  v  = 
is  the  radiation  wavenumber,  and  Jo  is  the  Bessel 
function  of  the  first  kind  of  order  zero.  The  terms 
ri  and  7*2  denote  the  distance  from  a  field  point 
to  a  source  and  to  its  image,  respectively,  and  are 
given  by: 

—  \j  (®p  ““  Xq)^  “1“  ( Up  ~  Vq  )2  ~{^p~-~Zq)^  (18) 

r2  =  sj (xp  -  ~Xq)2  +  (yp  -  j/,)2  +  (Zp  +  *,)2  (19) 

The  term  R  is  the  horizontal  distance  between  the 
field  and  source  points: 

R  =  \J (xp  -  ~Xq)2  +  (t/p  - 1/,)2  (20) 

Applying  the  hull  body  boundary  conditions  of 
Equations  13  and  14,  the  source  strengths  crga  and 
erg6  due  to  motions  of  ship-a  for  each  mode  k  are 
solved  by  satisfying  the  following  equations  simul¬ 
taneously: 

X  <rrto)ac^’8)ds(g) 

=  ni-f 77  on  Sa  (21) 

lOJe 

=  0  on  Sb  (22) 

Once  the  complex  source  densities  crga  and  erg6  are 
solved,  the  complex  radiation  velocity  potentials 
(f> g  can  be  determined  using  Equation  (16).  For 
radiation  due  to  motions  of  ship-b  while  ship-a  is 
stationary,  the  source  densities  crga  and  erg6  and 
potentials  can  be  solved  in  a  similar  manner. 

Added  Mass  and  Damping 

Once  the  radiation  potentials  are  known,  the 
added  mass  and  damping  terms  can  be  evaluated. 


The  radiation  force  component  j  acting  on  ship-a 
can  be  expressed  as: 

Ffa{x,y,z,t)  = 

£  [~Wk  -  C£A?£  -  <»/$  -  Cp$]  (23) 

k= 1 

where  /x|g  is  added  mass  due  to  motion  compo¬ 
nent  k  of  ship-  a,  is  damping  due  to  motion 
component  k  of  ship-a,  is  added  mass  due  to 
motion  component  k  of  ship-b,  and  is  damp¬ 
ing  due  to  motion  component  k  of  ship-b.  Based 
on  Equations  (8),  (10),  and  (23),  the  added  mass 
and  damping  coefficients  for  forces  on  ship-a  are: 


tfk  =  p{js 


The  added  mass  and  damping  coefficients  for  forces 
on  ship-b  are  evaluated  in  a  similar  manner. 


WAVE  EXCITATION  FORCE 

The  wave  excitation  force  on  a  ship-a  in  reg¬ 
ular  waves  can  be  expressed  in  terms  of  complex 
forces  due  to  incident  and  diffracted  waves: 

FjVa(t)  =  Re  [(//«  +  //>“)e-fa'«t]  (28) 
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where  fja  is  the  complex  incident  (Froude-Krylov) 
force  and  fPa  is  the  diffracted  wave  force. 


Incident  Wave  Force 

The  complex  incident  wave  potential  intro¬ 
duced  in  Equation  (2)  is  the  following  for  unit  am¬ 
plitude  waves: 


ship  hulls  for  the  complex  diffracted  potential  4>n 
are: 


d(j>p 

dn 

d(j)p 

dn 


<¥/ 

dn 

on  Sa 

(33) 

dfa 

dn 

on  Sb 

(34) 

The  diffracted  potential  is  represented  by  a  series 
of  sources  on  both  ships: 


=  JL  .  ekz+ik(x  cos /3+y  sin  (3)  £29) 

iuj 

where  g  is  acceleration  due  to  gravity  and  /?  is  the 
wave  direction  relative  to  the  ship  heading  (see  Fig¬ 
ure  1).  The  wavenumber  fc  of  the  incident  waves  is 
given  by: 

up 

k  =  —  (30) 

9 

The  relationship  between  the  wave  encounter  fre¬ 
quency  and  wave  frequency  is: 


<t>D(p)  =  /  0£>(g)<3(p;  q)dS{q) 

JS . 

+  /  c TbD{q)G(p;q)dS(q)  (35) 

Jsb 

where  crf^(q)  is  the  source  density  on  ship-a  and 
a bD(q)  is  the  source  density  on  ship-b.  The  Green 
function  G{jp\q)  of  equation  (17)  again  satisfies  all 
conditions  except  the  ship  hull  boundary  condi¬ 
tions.  Applying  the  body  surface  boundary  condi¬ 
tions,  the  diffracted  source  strengths  are  solved  to 
satisfy  the  following  equations  simultaneously: 


u)e  =  uj  —  kU  cos/3  (31) 


Substituting  the  complex  incident  wave  potential 
(j>j  into  equation  (8),  the  complex  incident  wave 
excitation  force  is: 


h 


la  _ 


—pcu  /  Im  [ <f>i ]  UjdS 
Jsa 

-p  f  Re  [V0/  *  V0*]  UjdS 
Jsa 

+i^ypuj  j  Re  [0/]  n^dS 

-p  J  Im  [V<£/  •  V4>„]  n?ds)  (32) 


+ls/^?£^rldS^  =  -d-st  -  onS“(36) 

+  at>D{q)r2^EiAdS{q)  =  on  Sb  (37) 

After  solving  the  diffracted  source  strengths, 
the  diffracted  potentials  are  then  computed  us¬ 
ing  equation  (35).  These  diffracted  potentials  can 
then  be  substituted  into  equation  (8)  to  get  the 
diffracted  wave  force  on  ship-a: 


The  above  equation  utilizes  the  relationship  be¬ 
tween  wave  frequency  and  wave  encounter  fre¬ 
quency  of  Equation  (31)  to  include  the  speed  de¬ 
pendent  terms  of  Equation  (8).  The  excitation 
force  on  ship-b  is  found  in  a  similar  manner. 


Diffracted  Wave  Force 

The  presence  of  the  two  ships  in  waves  will 
modify  the  incident  wave  field,  creating  a  diffracted 
wave  force  on  each  ship.  The  diffracted  wave  forces 
are  solved  using  a  source  distribution  method  sim¬ 
ilar  to  that  used  for  the  radiation  forces;  how¬ 
ever,  the  diffracted  wave  potential  only  needs  to 
be  solved  for  the  case  of  two  stationary  ships  in 
waves.  The  resulting  boundary  conditions  on  the 


fDa 


=  -p^e  J  Im  [<j>D]  UjdS 
~p  J  Re [V0z>  •  V<£s]  TijdS 

+iiypue  J  R e[<f>D]  rijdS 


L 


+pU  /  Im 


d<f>D 


dx 


-p  j  Im  [V<f>D  ■  V<£g]n“ds) 


(38) 


The  diffracted  forces  on  ship-b  are  solved  in  a  sim¬ 
ilar  manner. 
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COMPUTATIONS  FOR  TWO 
CIRCULAR  CYLINDERS 

The  computer  program  was  initially  validated 
using  results  from  References  3  and  4.  Figure  2 
presents  the  case  of  two  stationary  vertical  circu¬ 
lar  cylinders  floating  in  waves,  which  was  used  for 
comparison  with  computations  by  Kokkinowrachos 
[4].  Each  cylinder  has  radius  r  and  draft  equal 
to  1/2  r.  The  computations  presented  here  are 
for  a  separation  distance  between  vertical  axes  of 
5r.  The  wave  direction  is  along  the  x  axis,  which 
passes  through  both  cylinders.  For  the  present  dis¬ 
cussion,  the  cylinder  on  the  left  is  referred  to  as  ex¬ 
posed  cylinder-a,  while  the  cylinder  on  the  right  is 
referred  to  as  sheltered  cylinder-b.  Each  cylinder 
has  92  panels  for  computations. 


Cylinder-a  Cylinder-b 

(exposed)  (sheltered) 


Figure  2:  Two  Identical  Circular  Cylinders  Float¬ 
ing  in  Waves 

Figures  3  to  8  show  sample  computed  results 
for  the  two  cylinder  case  of  Figure  2  and  for  a 
single  cylinder.  Agreement  with  computations  by 
Kokkinowrachos  for  two  cylinders  is  very  good, 
suggesting  that  the  method  has  been  correctly  im¬ 
plemented.  Figures  3  and  4  indicate  that  the 
surge  excitation  force  on  the  exposed  cylinder  is 
influenced  significantly  by  the  sheltered  cylinder; 
however,  the  influence  on  the  sheltered  cylinder 
is  small.  For  heave  excitation,  interaction  effects 
have  only  a  small  influence  on  both  cylinders.  The 
heave  interaction  added  mass  and  damping  terms 
shown  in  Figures  7  and  8  appear  to  be  significant. 
Because  the  two  cylinders  have  identical  dimen¬ 
sions,  the  interaction  terms  and  A  are  equal 
to  /xfj  and  for  all  values  of  i  and  j. 

WAVE  EXCITATION  FORCES  ON 
TWO  SHIPS 

After  successful  validation,  the  program  was 
run  for  a  naval  Maritime  Coastal  Defence  Vessel 
(MCDV)  abeam  of  a  supply  vessel  (AOR),  with 


Dimensionless  Wavenumber  kr 


Figure  3:  Surge  Excitation  on  Exposed  Cylinder-a 


Dimensionless  Wavenumber  kr 


Figure  4:  Surge  Excitation  on  Sheltered 

Cylinder-b 


Dimensionless  Wavenumber  kr 


Figure  5:  Heave  Excitation  on  Exposed  Cylinder-a 
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Figure  6:  Heave  Excitation  on  Sheltered 

Cylinder-b 


Figure  7:  Heave  Added  Mass  Interaction  for 
Cylinders 


Figure  8:  Heave  Damping  Interaction  for  Cylinders 


both  ships  travelling  at  10  knots.  Table  1  gives 
dimensions  for  both  ships.  Using  the  notation  of 
Figure  1,  the  AOR  is  ship-a  and  the  smaller  MCDV 
is  ship-b.  For  the  present  computations,  midships 
of  the  MCDV  is  aligned  longitudinally  with  mid¬ 
ships  of  the  AOR  (i.e.  —  0).  Lateral  offsets  of 

25.44  m  and  500  m  for  the  MCDV  are  used,  giving 
a  narrow  gap  of  8.5  m  between  ships  for  the  first 
case  and  a  very  wide  separation  distance  for  the 
second  case. 

Each  ship  is  panelled  with  approximately  20 
panels  longitudinally  by  10  panels  girthwise,  re¬ 
sulting  in  200  panels  for  the  AOR  and  212  pan¬ 
els  for  the  MCDV.  The  maximum  nominal  panel 
dimension  is  approximately  8  m,  representing  the 
longitudinal  dimension  of  panels  on  the  AOR  hull. 
According  to  Sarpkaya  and  Isaacson  [5],  a  panel 
method  should  be  able  to  give  accurate  results 
for  wavelengths  greater  than  8  times  the  maxi¬ 
mum  panel  dimension;  thus,  the  current  panelling 
scheme  should  give  good  results  for  wavelengths 
greater  than  64  m,  representing  \/L &  >  1.2.  The 
panel  size  should  also  be  less  than  the  gap  between 
two  bodies  in  close  proximity,  which  is  satisfied  for 
the  current  cases. 


Table  1:  Dimensions  of  AOR  and  MCDV 


AOR 

MCDV 

Length,  L 

162.5  m 

52.0  m 

Beam,  B 

23.16  m 

10.76  m 

Draft,  T 

9.14  m 

3.4  m 

Block  coefficient,  Cb 

0.60 

0.56 

Beam  Seas 

Figures  9  to  16  give  sway  and  heave  excitation 
forces  acting  on  the  MCDV  and  AOR  from  the 
two  possible  beam  seas  directions.  Wave  length 
A  is  non-dimensionalized  by  the  MCDV  length  L b 
for  all  figures.  For  each  combination  of  ship  speed, 
wave  heading,  separation  configuration,  and  wave¬ 
length,  the  computation  time  was  approximately 
one  hour  on  a  workstation  with  speed  ratings  of 
85  MIPS  and  70  SPECmarks,  leading  to  a  compu¬ 
tation  time  of  approximately  60  hours  for  the  two 
headings,  two  separation  distances,  and  15  wave¬ 
lengths  presented  here. 

Figures  9  and  10  show  that  the  sway  excita¬ 
tion  on  the  MCDV  is  dramatically  reduced  when 
it  is  sheltered  by  the  AOR.  The  sheltering  effect  is 
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surprisingly  great  at  the  large  separation  distance 
of  500  m.  As  expected,  Figures  11  and  12  show 
that  the  MCDV  has  very  little  influence  on  the 
AOR  sway  excitation.  Figures  13  and  14  indicate 
that  interaction  effects  on  the  MCDV  are  signif¬ 
icant  for  heave  excitation  but  less  dramatic  than 
for  sway  excitation.  The  heave  excitation  on  the 
AOR  shown  in  Figures  15  and  16  varies  little  with 
proximity  to  the  MCDV. 


Dimensionless  Wavelength  A /L& 


Figure  9:  Sway  Excitation  on  MCDV  Abeam  of 
AOR,  =  90  degrees  (MCDV  Sheltered) 


Dimensionless  Wavelength  A/T& 


Figure  10:  Sway  Excitation  on  MCDV  Abeam  of 
AOR,  j3  =  270  degrees  (MCDV  Exposed) 


Head  Seas 

Because  ships  typically  head  into  waves  during 
replenishment-at-sea  operations,  additional  com¬ 
putations  were  performed  for  the  MCDV  abeam  of 
the  AOR  in  head  seas  at  a  speed  of  10  knots.  Fig¬ 
ures  17  to  20  show  heave  and  pitch  excitation  forces 
on  the  MCDV  and  AOR.  With  the  exception  of 
heave  excitation  on  the  MCDV  in  the  vicinity  of 
A  jL\y  —  2,  interaction  effects  are  very  small.  The 


Dimensionless  Wavelength  X/Li 


Figure  11:  Sway  Excitation  on  AOR  Abeam  of 
MCDV,  f3  —  90  degrees  (AOR  Exposed) 


Dimensionless  Wavelength  A /L& 


Figure  12:  Sway  Excitation  on  AOR  Abeam  of 
MCDV,  fi  —  270  degrees  (AOR  Sheltered) 


Dimensionless  Wavelength  A/L& 


Figure  13:  Heave  Excitation  on  MCDV  Abeam  of 
AOR,  (3  =  90  degrees  (MCDV  Sheltered) 
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Dimensionless  Wavelength  A/L& 


Figure  14:  Heave  Excitation  on  MCDV  Abeam  of 
AOR,  (3  =  270  degrees  (MCDV  Exposed) 


Figure  15:  Heave  Excitation  on  AOR  Abeam  of 
MCDV,  /?  =  90  degrees  (AOR  Exposed) 


Figure  16:  Heave  Excitation  on  AOR  Abeam  of 
MCDV,  /?  =  270  degrees  (AOR  Sheltered) 


large  increase  in  heave  excitation  for  the  MCDV 
could  be  due  to  either  wave  resonance  between  the 
two  vessels  or  to  numerical  problems  caused  by  the 
natrow  gap  between  the  vessels. 


r?Wb 

. fa _ 

pgaL^Bb 


Figure  17:  Heave  Excitation  on  MCDV  Abeam  of 
AOR,  Head  Seas 


Figure  18:  Heave  Excitation  on  AOR  Abeam  of 
MCDV,  Head  Seas 


CONCLUSIONS 

A  computer  program  has  been  developed  for 
computing  the  unsteady  hydrodynamic  forces  on 
two  ships  in  close  proximity.  The  aim  of  the  pro¬ 
gram  was  to  investigate  interaction  effects  during 
operations  such  as  replenishment-at-sea. 

The  mathematical  solution  satisfies  boundary 
conditions  on  both  ship  hulls  simultaneously  for 
the  radiation  and  diffraction  potentials.  The  lack 
of  lateral  symmetry  for  the  two  ship  case  causes 
the  number  of  panels  to  be  approximately  4  times 
greater  and  the  computational  effort  to  be  approx¬ 
imately  16  times  greater  than  for  a  single  ship. 
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Dimensionless  Wavelength  A/L& 


Figure  19:  Pitch  Excitation  on  MCDV  Abeam  of 
AOR,  Head  Seas 


Dimensionless  Wavelength  X/Lb 


Figure  20:  Pitch  Excitation  on  AOR  Abeam  of 
MCDV,  Head  Seas 


Comparisons  with  published  results  for  hydro- 
dynamic  forces  on  two  circular  cylinders  indicate 
that  the  computer  program  is  giving  correct  re¬ 
sults.  Computations  for  a  Maritime  Coastal  De¬ 
fence  Vessel  operating  abeam  a  supply  ship  show 
that  interaction  effects  are  very  significant  for  ex¬ 
citation  forces  in  beam  seas  but  are  small  in  head 
seas,  with  the  exception  of  a  narrow  wavelength 
range  for  which  a  resonance  phenomenon  appears 
to  occur. 

The  present  computer  program  has  poten¬ 
tial  to  assist  with  design  and  planning  for 
replenshment-at-sea;  however,  the  program  should 
first  be  validated  with  experimental  data  for  two 
ships  in  close  proximity.  Given  a  lack  of  experi¬ 
mental  data,  a  new  series  of  wave  tank  experiments 
would  likely  be  required  for  program  validation. 
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